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Abstract: Without demanding a specific form for the inflaton potential, we obtain 
an estimate of the contribution to the curvature perturbation generated during the 
linear era of the hybrid inflation waterfall. The spectrum of this contribution peaks 
at some wavenumber k = k*, and goes like k 3 for k <C h*, making it typically 
negligible on cosmological scales. The scale fc* can be outside the horizon at the end 
of inflation, in which case ( = — {g 2 — (g 2 )) with g gaussian. Taking this into account, 
the cosmological bound on the abundance of black holes is likely to be satisfied if 
the curvaton mass m much bigger than the Hubble parameter H, but is likely to be 
violated if m < H . Coming to the contribution to £ from the rest of the waterfall, we 
are led to consider the use of the 'end-of-inflation' formula, giving the contribution to 
£ generated during a sufficiently sharp transition from nearly-exponential inflation 
to non-inflation, and we state for the first time the criterion for the transition to be 
sufficiently sharp. Our formulas are applied to supersymmetric GUT inflation and 
to supernatural/running-mass inflation. 
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1. Hybrid inflation 

Hybrid inflation [ffl, |2], |3| ends with a phase transition known as the waterfall, which 
up to now has been studied only in special cases. This paper, which is a continuation 
of ||, provides a rather general treatment. We begin by defining the setup. 

1.1 Scales leaving the horizon 

An inflation model starts to make contact with observation only around the time 
that the observable universe leaves the horizon. The following description hybrid 
inflation is intended to apply to the subsequent era. 

Within the standard cosmology, the number iV bs of e-folds of inflation after the 
observable universe leaves the horizon satisfies [0] 

m -y^< N ^< i9 -yjq?L. (L1) 

(The time-dependence of H is ignored in this expression, which is usually a good 
approximation.) The upper bound corresponds to matter domination from the end 
of inflation to the epoch T = 1 MeV, with radiation domination thereafter until the 
observed matter dominated era, while the lower bound replaces the former era by 
one of radiation domination. 

The scales probed by observation of large scale structure (cosmological scales) 
leave the horizon during the first 15 or so e-folds after the observable universe. On 
these scales, the curvature perturbation ( is nearly gaussian with a a nearly scale- 
invariant spectrum V^(k) ~ (5 x 10~ 5 ) 2 . 

1.2 Hybrid inflation potential 

Our analysis applies to a wide class of hybrid inflation models. The essential features 
of the potential are captured by the following expression; 

V{<f>, X) = V + V(<f>) + \m\<\>) X 2 + \\X A (1-2) 
m 2 (0) ^gY-m^g 2 ^ 2 -^). (1.3) 

To have a perturbative quantum theory we demand g <C 1 and A 1. The inflaton </> 
is supposed to have zero vev and V(4>) is set to zero at the vev. We require V'((j)) > 
during inflation so that <fi moves towards its vev. The era of inflation with <ft < C is 
called the waterfall. 
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The requirements that V and dV/dx vanish in the vacuum determine Vo and 
the vev of the waterfall field \'- 

XL = m 2 /X, V = m 4 /4A. (1.4) 

It is necessary for our analysis that the waterfall field x h as the canonical kinetic 
term. For simplicity we will pretend that x is a single real field. At least within the 
Standard Scenario defined below, that cannot really be the case because it would lead 
to the formation of domain walls, located along surfaces at which x( x )^0 is trapped 
at the origin, which would be fatal to the cosmology. In reality, x will be replaced 
by a function of two or more real fields. So that there is only one effective degree of 
freedom in the x direction, we will demand that the function is invariant under some 
symmetry group of the action. Then the only change in our analysis for the realistic 
case would be the introduction of some numerical factors into the equations. In the 
realistic case the domain walls might be replaced by cosmic strings or monopoles, 
but in general the trapping of x will n °t occur and x( x , t) will everywhere approach 
its vev. 

The inflaton <p may also be replaced by a function of two or more real fields. 
If there is still only one effective degree of freedom the only change is again the 
introduction of numerical factors. In the opposite case of multi-field inflation, corre- 
sponding to a family of inflationary trajectories that are curved in field space, most 
of our analysis still applies if, by the onset of the waterfall, the family has collapsed 
to a single effective trajectory which has negligible curvature during the waterfall. 
To obtain powerful results we take the inflaton to have the canonical kinetic term, 
though much of our analysis would apply to, for instance, k-inflation 

Hybrid inflation was first discovered in in the context of single-field inflation 
[0, [|. It was given its name in ||, where the form ( L2]) was invoked for V(<f>,x) 



with V{4>) = rn^cf) 2 /2. With parameters chosen to give the Standard Scenario, and 
demanding also that 5(f) is responsible for the observed curvature perturbation, this 
gives spectral index n > 1 in contradiction with observation. Many forms of V(4>) 
have been proposed, which allow 5(fi to generate the curvature perturbation JF| [5| 
within the single-field inflation scenario. 

In our calculations we employ Eq. ( |1.2| ) for V(<f), x), without specifying the in- 
flaton potential V{4>). Minor variants of Eq. (|1.2| ) would make little difference. The 
interaction g 2 <p 2 x 2 might be replaced by 2 x 2+n /A n where A is a uv cutoff, or the 
term A% 4 might be replaced by x 4+n /A n - For our purpose, these variants are equiva- 
lent to allowing (respectively) g and A to be many orders of magnitude below unity. 

More drastic modifications of Eq. (|1.2| ) have been proposed, including inverted 
hybrid inflation || where is increasing during inflation, as well as mutated and 



smooth hybrid inflation || [T(J where the waterfall field varies during inflation. Also, 
the waterfall potential might have a local minimum at the origin so that the waterfall 
proceeds by bubble formation |TT], H|. Our analysis does not apply to those cases. 
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1.3 Standard Scenario 

By varying the parameters in the potential ( |1.2|) , one can have a wide range of 
scenarios that is still not fully explored. Most discussions of hybrid inflation make 
some assumptions, corresponding to what might be called the Standard Scenario. In 
this section we state those assumptions, which are made in the rest of the paper. 

Until x approaches its vev at the end of the waterfall, inflation is supposed to 
be nearly exponential (e# = \H\/H 2 <C 1) with Vq dominating the potential: 

3M^H 2 {t) = p(t) ~ V . (1.5) 

We take H to be constant during the waterfall, which is usually a good approxima- 
tion. Nearly exponential inflation requires 

<P 2 < 3M 2 H 2 . (1.6) 

Eqs. and give 

X\JM 2 ~ YlH 2 jm 2 . (1.7) 

It is usually supposed that x vcv <C Mp corresponding to m 3> H. (In particular, 
GUT inflation jTJ, |TJ] takes x to be a GUT Higgs field with Xvcv ~ 10" 2 M P .) One 
sometimes considers Xvev roughly of order Mp corresponding to m roughly of order 
H (supernatural |14[ and running mass ]15] inflation). There do not seem to be any 
papers considering x vev 3> M P which would correspond to m ^ H. 
Using Eqs. ( |1.4|) and ( [L.5| ), the requirement A 1 is equivalent to 



m/H < ^Mp/H (A<1). (1.8) 
Successful BBN and the upper bound on the tensor perturbation require 

10" 42 < H/M P < 1(T 5 (BBN and tensor). (1.9) 



The upper part of the range is favoured, especially [[RJ because we deal with hybrid 
inflation. 

One usually requires Mp but we will just invoke the weaker requirement* 1 

<p c = m/g < M P . (1.10) 

If <p is big enough we have m 2 (0) > H 2 . Then we assume that x vanishes up 
to a vacuum fluctuation which is set to zero. If 4> is small enough, m 2 (0) < —H 2 . 
Hence there is a 'transition' regime with |m 2 (0)| H 2 . If the transition takes 
several Hubble times, the quantum fluctuation of x w ih be converted to a classical 
perturbation, with spectrum ~ (H/2n) 2 on all scales leaving the horizon during the 



* 1 This is also invoked in our earlier paper jl but note that Eq. (80) of B has a typo. 
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transition. To avoid this the transition should take less than a Hubble time or so 
(fast transition). 

The waterfall starts at m 2 (<f)) = which is in the middle of the transition. During 
the waterfall the vacuum fluctuation of \ is converted to a classical field x( x ; t), with 
X 2 moving towards Xvev The waterfall ends when x 2 (x, t) ~ Xvev an d inflation is 
supposed to end then because V((p) is not supposed to support inflation without the 
additional term Vq. 

Regarding <ft, we require that it decreases monotonically before the waterfall, and 
afterward for as long as it affects the evolution of \. This assumption is not at all 
trivial, because V(<p) may steepen as (f> decreases, causing (f> to oscillate about the 
origin. The evolution of x has yet to be studied for that case, which occurs in part of 
the parameter space for some well-motivated forms of the potential, including GUT 
inflation [|l^, [H| and running- mass inflation f[5 |. 



1.4 The waterfall 

During the waterfall we need to consider both (f> and \. Taking both fields to live 
in unperturbed spacetime (ie. ignoring back-reaction) the evolution equations during 
the waterfall are 

+ 3 #0 + _ V 2 = -sY0 (1.11) 
X + SH X + m 2 (<P) X - V\ = -Ax 3 - (1-12) 

We assume that the waterfall starts with an era during which Eq. (|1.12|) can be 
replaced by* 2 

Xk + 3#Xk + [{k/af + m 2 {t))] Xk = 0, (1.13) 

with m 2 (t) is independent of x- We call this the linear era, and it will be our main 
focus. Regarding Xk(^) as an operator J|, ||, its mode function Xkif) also satisfies 
Eq. (|1.13|) . We will see how Xk{t) grows exponentially for suitably small k, generating 
a classical quantity Xk(t)- Keeping only the classical modes, we arrive at a classical 
field x- 

During at least the first part of the linear era, m 2 (0) depends significantly on 
(p. Then the right hand side of Eq. (|1.11|) has to be negligible so that m 2 (t) can 
be independent of x- With that condition in place we just have to worry about the 
perturbation 5(f) that is generated from the vacuum fluctuation. 

If the linear era of the waterfall takes no more than a Hubble time or so, 5(f) 
can be completely eliminated by taking the spacetime slicing to be one of uniform 
(f). But if the waterfall takes much more than a Hubble time, new contributions to 
5(f) are generated as each scale leaves the horizon. To avoid this quantum effect on 



# 2 This assumption implies some lower bound on |0| but it is not clear how to calculate the bound 

a hi. 



- 5 - 



the evolution Xi we have to assume that the new contributions to 5(f) are negligible. 
Then we can again choose the slicing of uniform 0.# 3 

For the threading of spacetime, we choose the comoving worldlines (those moving 
with the fluid, so that a fluid element has zero momentum density). The gradients 
of both and (as we shall see) the classical field x are small compared with their 
time derivatives. If they vanished, the comoving worldlines would be free-falling and 
orthogonal to the slicing, and we could choose the time coordinate labeling the slicing 
to be proper time along each thread. We assume that the gradients are small enough 
to make that choice possible to an acceptable approximation. This completes the 
definition of the gauge in which the classical field x( x > t) is defined. 

We also need to justify the use of Eq. ( |1.13|) for the mode function Xhi before 
the classical quantity Xkif) is generated. As we will see, Eq. ( |1.13| ) is needed for 
that purpose only for modes that are well inside the horizon during this time and (at 
least with the approximation of Section |5]) only for much less than a Hubble time. 
That being the case, we can ignore the second term of Eq. ( ET3 ) and set a equal to a 
constant so that Eq. ( |1.13| ) becomes a flat spacetime equation in which back-reaction 
is negligible. 



2. Waterfall field during the linear era 



2.1 Evolution of x 

With H constant, we can use conformal time rj 

d 2 (axk) 



with 



drj 2 



-1/aH to write ( |1.13| ) as 



+ ^k a Xv = 0, 



(2.1) 



ul{rj) = r + a 2 m 2 (t) 



m 2 



m 2 {t) -2H 2 



m 2 (t) = g 2 (j) 2 (t) -m\ (2.2) 



For sufficiently small k, we can set uj\ 



J l=o 



a 2 m 2 . 



Then uj\ switches from 
positive to negative before <fi = <p c (but much less than a Hubble time before, by 
virtue of our fast transition assumption). For k 2 > the switch is later. For the 
scales that we need to consider, we assume that there are eras both before and after 



the switch when uj\ satisfies the adiabaticity condition d\uj k \jdr] <C \uj 



21 
fcl- 



# 3 Since we neglect the new contributions to Scj) during the linear era (while m 2 (4>) depends 
significantly on 0), we neglect also their effect on the spectrum V^(k). That is presumably a 



good approximation if Vq^ given by Eq. (5.12) is much less than the contribution V( nil (k) that we 
are going to calculate. That will probably be the case if V( nD (k) is big enough to form black holes, 
which is our main concern. It will not be the case if we deal with one of those exceptional inflation 
models where V^Ak), on the scales leaving the horizon during the linear era of the waterfall, is 



itself big enough to form black holes. Then we have the opposite situation: Eq. (5.12) will be valid 
if V^ nil is not big enough to form black holes. 
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Taking Xk to be an operator, its mode function Xk satisfies Eq. (|2.1|) . During 
the adiabaticity era before the switch we take the mode function to be 



axk = [2uk(v)) 1/2 exp ( -i / u k {r))dr) 



which defines the vacuum state. During the adiabaticity era after the switch 



axk ^ (2\u h (rf)\) 1/2 exp ( / \u k {r})\dq 

l m (k) 



(2.3) 



(2.4) 



where the subscript 1 denotes the beginning of the adiabatic era. The displayed 
prefactor is exact M only if m 2 {t) oc t and H(t -fjcl. As we are about to see, 
Xk grows during this era and we call it the growth era. 

During the growth era, the adiabaticity condition is equivalent to the two con- 
ditions 



2H 

\m(t)\ 

d\m(t)\ 



\m{t)\ 2 dt 



< 



k 



1 



a(t)\m(t)\ 



a(t)\m(t)\ 



3/2 



3/2 



(2.5) 
(2.6) 



The growth era begins when both conditions are first satisfied. 

The first condition implies 3> H so that — |w(t)|, and it can hold 

only if m ^> H. For k ^ a(t)\m(t)\ we have \uk\ — |wfc=o| — a(t)\m(t)\. With 
Eq. (|2.4|) this gives Xk — |^fc|Xfc/°- At = Eq. (|2.4|) becomes 



(2.7) 



Ignoring the relatively slow time-dependence of the prefactor, we have as a rough 
approximation 

Xk=0 {t) ~ (2a?|m(t 1 )|)~ 1 / 2 exp ( [ dt\m{t)\] . (2.8) 



Xk =o(t) ~ (2a 3 |mWI)" 1/2 exp / dt\m(t)\) . 



In the regime k <C a|m(t)| we have 



\tut\ — a\m{t)\ I 1 



giving 



where 



2a 2 \m 2 (t)\ 

Xk(t)^Xk=o(t)e- k2 / 2k ^\ 
dt 



l a 2 \m{t)\ 



(2.9) 

(2.10) 
(2.11) 
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By virtue of Eq. (|2.6|) , the change in in time |m(t)| 1 is small and so is 

the change in a. Defining 

t 8tart =ti + |m(ti)|-\ (2.12) 

we get 

MWt)^a(*i)M*i)l- (2.13) 
To avoid the divergence in k*{t) at t — t\, we will regard t star t as the start of the 
growth era rather than t\. After t s tart, K(t) decreases while a\m(t)\ increases. We 
assume that k 2 (t) <C a 2 \m 2 (t) \ , except for a brief era near the beginning of the growth 
era that can be ignored.* 4 Then Xk(t) at fixed t falls exponentially in the regime 
k*{t) < k < a\m(t)\ and significant modes have k < k*(t). 

The number of e-folds of growth is N(t) = H{t — t stait ). We denote the end of 
the linear era by a subscript 'end'. If N(t en d) < 1, k*(t) falls continuously. If instead 
N(t en d) ^> 1, the exponential increase of a causes k*(t) to level off after N(t) ~ 1. 
Using H |^(^)| < m , we learn that in any case 

K(t) \ 2 m 



1 « ,. y \rr < 77 • ( 2 - 14 ) 

This tells us that the scale k*(t) is shorter than the scale leaving the horizon at the 
beginning of the waterfall. Since we assume that cosmological scales are outside the 
horizon at this stage, is shorter than any cosmological scale. Dividing both 

sides by exp(2iV(t en d)) we have 

^ , <<(^)^(I) 1/2 ^- < 2 - 15 > 

This tells us that scale k*(i) (and in particular, its final value fc*(t e nd)) can be far 
outside the horizon at the end of the linear era. 

2.2 Classical field x(x,t) 

During the growth era the mode function Xk nas constant phase (zero with our 
convention), which means that Xk(^) oc Xk{t) can be be regarded as a classical field. 
The significant modes have k < fc*(t) a(t)\m(t)\. Because when for each mode. 
This means that the continuous creation of new classical modes, occuring for each 
mode when uj\ becomes negative at k ~ a(t)m(t), can be ignored. 

For the significant modes, Xk/Xk — a(t)\m(t)\. The classical field has approxi- 
mately the same behaviour;* 5 

X{x,t)~\m{t)\ x {x,t). (2.16) 



# 4 This assumption holds within the approximation of Section |[ 

# 5 This behaviour breaks down near any locations with x(x, t) = 0. To discuss them we would 



have to extend the discussion to a multi-component \ as mentioned at the end of Section 1.2. We 
assume that if they exist, they are rare enough to be ignored. 
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Since <C o|m(t)|, the gradient of x is small compared with its time derivative. 
The spectrum of x is 

V x (k,t) EE (fc 3 /27T 2 )P x (fc,t) = (fc 3 /27T 2 )| Xfc (t)| 2 . (2.17) 

Using Eq. ( [2~ToD the mean square (spatial average) of x 2 is jD 

(X 2 (t)) = J fv x (k,t) = 7^)3 / *kP x (k,t) = (2n)-^P x (0,t)kUt), (2.18) 

where P x (0,t) = \ X k=o(t)\ 2 is given by Eq. Q. 
We denote the perturbation in x 2 by ^X 2; 

5 X 2 (x,i)= X 2 (x,t)-( X 2 (t)>. (2.19) 

The convolution theorem gives |1| for P^ = (2ir 2 / k 3 )Vs x 2 

P Sx <k,t) = 7^53/ ^'P x (A; / ,t)P x (|k-k'|,t). (2.20) 

For fc < this gives @] 

n x2 (A:,t) = ^(x 2 W) 2 [A:A,(t)] 3 . (2.21) 

For ^> it gives 

V Sx 2(k,t) = 2( X 2 (t))V x (k,t), (2.22) 
which falls exponentially at fixed t. 

2.3 End of the linear era 

At each location, the linear equation fl2.16p ceases to be valid around the time when 



X 2 (x, t) achieves some value Xn\- This time is given by 

X 2 (x,t nl (x)) = X 2 1 . (2.23) 

We will take the linear epoch to end at a time t en d, such that the fraction y > of 
space with x 2 (x, t) > xL\ is small. We will see that the probability distribution of 
x(x, t) is gaussian, and using the approximation erfc(x) ~ exp(— x 2 ) we have 

X 2 1 ^ln(l/y > )(x 2 (tend)}- (2.24) 

According to this equation t end is not very sensitive to the choice of y > , and for 
estimates we will take ln(l/y>) ~ 1. 

If the linear era lasts for long enough, we will have m 2 (t en d) — — m 2 (ie. 0(t en d) "C 
4> c = m/g). In that case the right hand side of Eq. (|1.11|) is irrelevant, and the linear 
era ends only when the right hand side of Eq. ( |1.12[ ) becomes significant. This gives 

Xl\ ^ XL = 12M 2 # 2 /m 2 . (2.25) 
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Now suppose instead that |ra(t e nd)| m (ie. 0(t en d) — 4>c)- Then the linear era 
will end when the right hand side of Eq. ( |1.11|) becomes significant, provided that 
the right hand side of Eq. ( pL .12[ ) is then still negligible which we are about to show 
will be the case. If is still slowly rolling so that 3H(p(t) = —V, the right hand side 
of Eq. ( |1.11| ) becomes significant, when x 2 ~ 3i/|0(O)|/gm. But it may be that g 2 x 2 
has become of order if 2 first, causing to oscillate about the origin. Including both 
possibilities we have* 6 



/3ff0(O) H 



x l~ min J^lZ±l± 7 ^y (2.26) 
y gm g 2 



It follows from Eqs. (|1.10|) and (|1.6| ) that the right hand side of Eq. ( |2.26|) is much 



less than Xvev> ma king the right hand side of Eq. ( P-.12j ) insignificant as advertised. 
2.4 Energy density and pressure of x 

We have seen that the gradient of x is negligible compared with its time-derivative, 
and in our adopted gauge the gradient of vanishes. Ignoring the gradient, the 
energy density and pressure are p = p^ + p x and p = p^ + p x where 

P^t) = l y o + V<M + (2.27) 
P4>V) = -[V +V(<f>)] + ~<j> 2 , (2.28) 

Px^-\\ m \ t )\x 2 + \x 2 (2-29) 
~ (2.30) 

Px^ \W(t)\x 2 + \x 2 (2.31) 

~ \m 2 (t)\ X 2 ^X 2 - (2.32) 

In an unperturbed universe the energy continuity equation holds; 

p{t) = -ZH(t)(p(t)+p(t)). (2.33) 

To the extent that spatial gradients are negligible it holds at each location. With a 
generic choice of the slicing, denoted by a subscript g, we have 

p g (x, t) ~ - 3 ^M (pg(X; t ) + Pg (x, t)) , (2.34) 

where a g (x, t) is the locally defined scale factor. We are working in the gauge defined 
in Section [D|, which means that t is proper time and we can choose a(x, t) = a(t), 
the unperturbed scale factor. We therefore have 

p(x, t) ~ -3H(t) (p(x, t) + p(x, t)) . (2.35) 



^ 6 With the potential V((f>) = ^m^cj) 2 the right hand side is min{m^, H 2 } = m| which means that 
only the first possibility exists. The existence of the second possibility for a more general potential 
was missed in 101. 
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We are dealing with the linear era, which means that the right hand sides of 
Eqs. Ql.lip and ( |1.12|) are negligible. With its right hand side negligible, Eq. ( |1.11| ) 



describes a free field which means that it satisfies the energy continuity equation by 
itself; 

P4> ~ -3H ( P<t> + p<j> ) = -3H0 2 . (2.36) 
The same must therefore be true for \] 

p x (x, t) ~ -3if (p x (x, t) + p x (x, t)) = -3H X 2 (x, t). (2.37) 

At each location we have* 7 

X(x, t) + 3#x(x, t) + m 2 (t) X (x, t) ~ 0. (2.38) 



Using this equation to differentiate Eq. ( 2.29 ) we find 



p x (x, t) ~ -3ff(p x (x, f) + p x (x, t)) - ^^^X 2 (x, «)■ (2-39) 

The second term of the right hand side violates the energy continuity equation. This 
apparent inconsistency between the field equations and the energy continuity equa- 
tion occurs because the effect of the interaction term g 2 <f) 2 x 2 is dropped in Eq. ( |1 . 1 1| ) 
(ie. the right hand side is set to zero) but kept in Eq. Q1.12 ). 



If we demand approximate consistency between the field equations and the en- 
ergy continuity equation, we need the second term of Eq. ( |2.39| ) to be much smaller 
than the first. That condition is equivalent to 

^ « \m(t)\H, (2.40) 



which is stronger than the adiabaticity conditions (2.5) and ( |2.6|) (with k a\m(t)\). 



But there is no need to impose this stronger condition, because the near cancellation 
of the two terms of p x makes it unreasonable to expect the approximate evolution 
( |2.38 ) of x to give even an approximate estimate of p x . By contrast, the right hand 



side of the energy continuity equation has no cancellation so that it can be used to 



evaluate p x . Invoking Eq. (|2.16|) , we find 

p x (x,t)^-^|m(t)| X 2 (x,t). (2.41) 
2.5 Justifying the neglect of \ 

It has been essential for our discussion that the spatial average of x(x, t) is negligible. 
That is the case in a sufficiently large volume, because x( x 5 1) is constructed entirely 



# 7 In the present context Eq. (2.38) can be replaced by Eq. ( |2.16| ) by virtue of the adiabaticity 



condition. But in Section ^ we drop the adiabaticity condition. 
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from the Fourier modes. But to make contact with cosmological observations we 
should consider a finite box, whose (coordinate) size L is not too many orders of 
magnitude bigger than the size of the presently observable universe ||19|| . Denote the 
average within the box by x we have 

X (x)~x + X>(x), (2.42) 

where the Fourier modes of x> satisfy kL > 1 so that 

(Xl) = jf %V x {k). (2.43) 

The average within the box comes from modes with kL < 1, and for a random 
location of the box the expectation value of x 2 is 

f L 1 dk 

(X 2 ) ~ J -V x {k). (2.44) 

To justify the neglect of x we need (x 2 ) -C (x>)- I n our scenario, where V x (k) 
peaks at a value k*(t), this is equivalent to Lk*(t) 3> 1. That is satisfied because the 
scale is supposed to be much smaller than the observable universe. To have 
(x 2 ) > (x>) we wou ld presumably have to allow the transition from m 2 (t) = —H 2 to 
m 2 (t) = H 2 to take at least several e- folds so that it can generate a contribution to % 
that has a flat spectrum. Then, if the flat spectrum generated during the transition 
dominates, one would have 

(X 2 )/(X>) - AWe/iW, (2.45) 

where iVbefore (iV" a f t e r ) is the number of e-folds of transition before (after) the observ- 
able universe leaves the horizon. 



3. Contribution ( x to the curvature perturbation 

We write the contribution to the curvature perturbation that is generated during the 
waterfall as ( x = Cim + Cnb where the first term is generated during the linear era, and 
the second is generated afterward up to some epoch just after inflation has ended. 

The curvature perturbation is C( x )^) = 51na(x, t), where a(x, t) is the locally 
defined scale factor on the spacetime slicing of uniform p. As in Eq. (|2.35| ), the 
spatial gradient is supposed to be negligible, which in general requires smoothing on 
a super-horizon scale. Using that equation, we see that the change in ( between any 
two times t 1 and t 2 is 

C(x, t 2 ) - C(x, t x ) = 5N(x, t a , t 2 ), (3.1) 
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where iV(x, ^1,^2) where iV is the number of e-folds between slices of uniform p. 
Working to first order in 5p,# 8 we will use this result to calculate (j in , and then see 
how it might be used to calculate Cni- 

We note in passing that an equivalent procedure is to integrate the expression 

Hit) 

C(x,t) = ~ (t) )p (t) foWx,t), (3.2) 

where 

pit) 

5p n ad(x, t) = 8p(x, t) - -—5p{bfx, t). (3.3) 

pit) 

3.1 The contribution (j in 

During the linear era, the gradient of p is negligible without any smoothing. We are 
working in a gauge where 5(f) = so that p(x, t) = p x (x,t) + p^t). Ignoring the 
inhomogeneity of the locally defined Hubble parameter, 

Clin(x, t)=H [5t(x, tend) " <*t(x, £ sta rt)] , (3.4) 

and 

M(x, t) = - H 5p ^ t] = i . (3.5) 
p(t) 3(^(t)) + 2 (t) 



Using Eq. (PUD , 



™ (X ' t) -3(x 2 (t))(^(t)) + ^(t) (X 2 W) ^ 



= 1 g (X 2 W) *X 2 (M) 

2|m(O|(x 2 (*)) + 2 (t) (* 2 (*)> 

Using Eq. (p.21|) , we have for <C 



(3.7) 



5H ; v2Kt)i; V(x 2 (t)) + 2 (t); vm*)/ 1 ; 

We assume that |5t(t e nd)| ^ |^(t start )|, which will be justified within the approx- 
imation of Section [|. Then we have 

Cn n (x) = -H5t(x,t end ) (3.9) 

H (X 2 (tcnd)) 5x 2 (x,tend) 



2|m(W)| (X 2 (tend)) + 2 (tend) (% 2 (W)) 



(3.10) 



# 8 A second-order calculation of ( is needed only to treat very small non-gaussianity corresponding 
to reduced bispectrum |/nl| ^ 1- On cosmological scales, such non-gaussianity will eventually be 
measurable (and is expected if £ comes from a curvaton-type mechanism || ) . But there is no hope 
of detecting such non-gaussianity on much smaller scales. 
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The inhomogeneity of H is indeed negligible because it generates a contribution 



6H(t,x)dt = — r" 6p(x,t)dt ~ (3.11) 



'Wt 2 ftiw ' 2p|m(t end )| 

This is much less than H6t(t en d) in magnitude, because \m(t)\ 3> H and |p| <C if p. 
Using Eq. (|3.8|), we have for k <C k*(t end ) 

AFP (k)=( - V7 {ntcnd)) ) 2 ( k W( k 
Clinl V2|m(W)|>/ \(x 2 (t end )) + ^(t cnd )J \h(t end )J ^\h(t cnd ) 

(3.12) 

In the opposite regime k 3> k*(t end ), Cim is negligible because 5\ 2 is. Therefore, P^ lin 
peaks at k*(t ead ) with the value 

V2|m(t cnd )i; V(x 2 (t C nd)) + 2 (w)y 

If m 2 (t en d) — — m 2 , Eq. (|2.25|) holds and we have 

(X 2 (t cnd ))/4> 2 (t cnd ) ~ 12M|#70 2 (t cnd ) > 1. (3.14) 
(The inequality follows from Eq. fll.fip ). Then Eq. ( |3.12| ) simplifies to 

k 



^<*> = © ( 



k* (t 



end , 



(3.15) 



3.2 The contribution £ 



nl 



Let us estimate the number of e-folds N n \ after the end of the linear era. At tend; X 
is increasing exponentially. Soon afterward it starts to affect 0, driving it towards 
zero. We therefore expect m 2 (t) to quickly approach — m 2 after t end (if it is not there 
already) , restoring at least approximately the linear evolution of x 2 ■ Then Eq. ( |2.16|) 
will hold with ~ m giving 

iV n i ~ (H/m) ln(Xvev/Xend)- (3.16) 

If the linear era ends only when the right hand side of Eq. ( [L.12 ) becomes important 
we have hi(x V cv/Xcnd) ~ 1, giving N n \ ~ H/m <C 1. But if it ends when the right 
hand side of Eq. flLllQ becomes important we may have ln(x V ev/Xcnd) 1 which 
allows iV n i > 1. 

Now we consider the contribution £ nl , that is generated between t en d and some 
time t<i just after inflation has everywhere ended. To calculate it we need to smooth 
on a super-horizon scale. Then we can use the 5N formula which gives 

Cni(x,t) =H6t 12 {x), (3.17) 
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where the initial and final slices both have uniform p and £12 (x) is the proper time 
interval between them. 

At each location, the linear era ends at the epoch £ n i( x ) given by Eq. (|2.23| ). 
At this epoch there is nearly-exponential inflation, and inflation ends at some later 
time tnoninf(x). If At(x) = t n oninf(x) — t n i(x) is sufficiently small it can be taken to 
correspond to a spacetime slice of negligible thickness. Then 5ti 2 is given by the 



'end-of- inflation' formula 2C 



<ft 12 (x) ~ (3.18) 

Pi tend J 

where Sp(x) is defined on the slice. The addition of Cni to Cim corresponds to taking 
the final slice of the SN formula be the transition slice, instead of a slice of uniform 
P- 

This equation is valid to first order in 5p. To derive it we take the separation 
between the initial and final slice to be not much bigger than is needed for them to 
enclose the transition slice. Then we can take the unperturbed quantity p(t) to have 
a constant value both during inflation and non-inflation. This gives 

5t 12 (x) ~ 5p(x) (— - . (3.19) 

\ Pinf Pnoninf / 

Since |pi n f| is evaluated during nearly-exponential inflation, it is much smaller than 
I Pnoninf I leading to Eq. ( |3.18| ). 



We are defining 5p x on a slice of uniform p^ and 5p is defined on a slice of uniform 
X- The time displacement from the first slice to the second slice is is —5p x /p x , which 
means that 5p = 5 p x P0 (£ e nd) / Px(£end) . Putting this into Eq. ( ft.lSQ we get 

Cnl (x) /Clin (x) ~ p^(tend)/p X (£end) ^ ^ (tend) / jf ( W) • (3.20) 

We are only interested in the case that this ratio is > 1. Then the inclusion of Cni 
corresponds to omitting the middle term of Eq. ( 3.10Q . From Eq. ( 3.14 ), this case 
can occur only if |m(t en d)| ^ m - 



Now comes a crucial point. From the derivation of Eq. ([3.18 ), it is clear that the 



criterion for its validity is At(x) |<5ti2(x)|, at a typical location. (In words, the 
thickness of the transition slice is negligible compared compared with its warping.) 
This simple remark has not been made before, and consequently it has not been 
checked whether the criterion is satisfied. 

In our case, HAt(x) is given by Eq. ( |3.16|) with xlnd replaced by Xni( x )- Since 



that quantity appears only in the log the change will not have much effect, and we 
will have i?At(x) ~ iV n i at a typical location. On the other hand, the typical value of 

1 /2 

|Cni(x)| = -ff |5ti2(x)| is (k) where k < aH is the smoothing scale used to define 
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Cni- The criterion for Eq. ( |3.18 ) to be valid is therefore N a \ <C V^(k). In the regime 



of interest 2 (t en d)/x 2 (^end) 3> 1, this criterion becomes 

|m(W)l In (^) « ( -^-) ^ .. (3.21) 

m \XendJ \k*{t end )J 

Whenever the criterion ( |3.21| ) is not satisfied, the calculation of Cend that we have 
described does not apply. 

3.3 Other uses of the 'end of inflation' formula 

Our use of Eq. ( |3.18| ) to evaluate ( n \ is quite different from its usual applications 



PH , p2] , |23"1 , p4}] . In those applications, the field causing <5p(x) has a nearly flat 
spectrum, leading to a nearly flat V^ ul (k) that can give a significant (even dominant) 
contribution to V^(k) on cosmological scales. Since V^ 2 {k) ~ 5 x 10~ 5 on these 



scales, Eq. ( |3.21| ) on cosmological scales becomes 

iV tra „ < V X £ < 5 x 10" 5 , . (3.22) 

where iV tran now refers to the duration of the transition slice in the scenario under 
consideration, and C12 = H5ti 2 is the contribution to (. 

Most of the other applications |20, 21, 22| consider hybrid inflation, with the 



transition slice the entire hybrid inflation waterfall. Of course their setup is different 
from ours because they introduce a third field, the one that generates 5p in Eq. (|3.18|) . 
In these cases, iV tran in Eq. ( |3.22|) becomes the total duration of the waterfall. We see 
from Eq. ( |3.16| ) that it cannot be much less than H/m, which means that Eq. ( |3.22j ) 
needs H/m < 5x 10~ 5 . Since we need (H/m) 2 3> H/Mp (corresponding to A <C 1), 
this requires a fairly low inflation scale H 10 _9 Mp. 

An alternative possibility [^3j is for the transition slice to be at the end of 
thermal inflation [^, E6[ ^7\, |2B[. (Thermal inflation is is a few e- folds of inflation 
occurring typically long after the usual inflation, which is ended by a thermal phase 
transition.) Then we expect roughly AiV ~ H/m, where m is the tachyonic mass of 
the field causing the end of thermal inflation. This criterion ( |3.22|) is satisfied by the 
usual realizations of thermal inflation. Further possibilities for the transition slice 



are considered in 24 . 



4. Effect of C x 

4.1 Cosmological black hole bound on V$ 

The most dramatic effect of ( would be the formation of black holes. This places 
an upper bound on which we now discuss taking on board for the first time the 
non-gaussianity of 
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The bound that we are going to consider rests on the validity of the following 
statement: if, at any epoch after inflation, there are roughly spherical and horizon- 
sized regions with ( significantly bigger than 1, a significant fraction of them will 
collapse to form roughly horizon-sized black holes.* 9 The validity is suggested by 
the following argument: the overdensity at horizon entry is Sp/p ~ (, and if it is of 
order 1 then 5p ~ p = 3MpH 2 . The excess energy within the Hubble distance H^ 1 
is then M ~ H~ 3 p ~ Mp/H, which means that the Hubble distance corresponds 
roughly to the Schwarzchild radius of a black hole with mass M. The validity is 
confirmed by detailed calculation using several different approaches, as summarized 



for instance in [2S 



Before continuing we mention the following caveat. Practically all of the liter- 
ature, as well as the simple argument just given, assumes that ( within the region 
is not very much bigger than 1. Then the spatial geometry within the region is not 
too strongly distorted and the size of the black hole is indeed roughly that of the 
horizon. In the opposite case, the background geometry is strongly distorted and the 
wavenumber k defined in the background no longer specifies the physical size of the 
region at the epoch aH = k of horizon entry [30||. An entirely different discussion 



would then be necessary, which has not been given in the literature. As the opposite 
case does not arise in typical early- universe scenarios we ignore it. 

We are interested in the case that V^(k) has a peak at some value fc pe ak, and we 
assume that the width of the peak in In k is roughly of order 1 so that 



POO 

(C 2 ) = / V ( (k)dk/k ~ V ( (k pe , k ). (4.1) 
Jo 



Regions with £ > 1 that might form black holes will be rare if P^fcpeak) is not too 
big. Observation demands that the regions must indeed be rare, because it places a 
strong upper bound on the fraction of of space that can collapse to form horizon-sized 
black holes, on the assumption that the collapse takes place at a single epoch as is 
the case in our scenario. A recent investigation of the bound is given in |29|| , with 



extensive references to the literature. The bound depends on the epoch of collapse. 
Denoting it by (3 it lies in the range 

10~ 20 < P < 10- 5 . (4.2) 

To bound P^fcpeak), we shall require y < (3 where y is the fraction of space with 
( > (c and ( c is roughly of order 1. 

The fraction y can be calculated from (( 2 ) if we know the probability distribution 
of C( x )- The standard assumption is that it is gaussian. Then 

y= 1 -erk(( c / v f2(e)), (4.3) 



# 9 We are choosing the background scale factor a{t) so that the perturbation ( = 5(lna(x, t)) has 
zero spatial average. 



-17- 



and using the large- a; approximation erfc (x) ~ e 3=2 / y/nx ~ e a ' 2 we find 

n(^Pca k )^(C 2 )<Cc 2 /21n(l//). (4.4) 



For the range Q this gives (with Cc ^ 1) ^ c (A; pcak ) < 0.01 to 0.04. 

But Clin given by Eqs. fl2.19|) and ( |3.10| ) is actually non-gaussian, of the form 



C = -{g 2 - (g 2 )). 



(4.5) 



With this form, there is no region of space where ( > (g 2 ), and 1/ < 1 now im- 
plies some bound (g 2 ) — Cc ^ Cc which is practically equivalent to (g 2 ) < Cc- This 
corresponds to V(;{k pcak ) < 2 

* (C 2 ) = 2(/) 2 < 2C 2 < 2. (4.6) 

For completeness, we see what happens if ( = + (g 2 — (g 2 )) with g gaussian. 
(This might be the case [^TJ if C is generated after inflation by a curvaton-type 
mechanism.) The we have 



^c(W) ~ (C 2 ) = 2(^ 2 ) 2 < 2 



Cc 



21n(l/y). 



(4.7) 



which gives P c (fc peak ) < 2 x 10" 4 to 2 x 10~ 3 . 

In all three cases, the bound on ^(fcpeak) is very insensitive to / which means 
that it depends only weakly on the value of /3. Turning that around though, the 
black hole abundance is very sensitive to Pf(/c pea k) which suggests that fine-tuning 
of parameters will be needed to get an eventually observable (yet presently allowed) 
abundance. 

If the peak has width A In k different from 1, (C 2 ) — "Pf (A; pea k) A In k. If A In k <C 
1 this weakens the bound on P^(/c pea k) by a factor (A In A;) -1 , but such a narrow peak 
is not generated in typical scenarios. If instead A In A; 3> 1, one might think that the 
bound on P^(fc pea k) is strengthened by a factor (Alnfc) -1 , but that conclusion is too 
hasty because the observational bound ( |4.2|) refers to the formation of horizon sized 
black holes at a more or less definite epoch whereas the broad peak will lead to the 
formation of such black holes over A In k Hubble times. The value of (C 2 ) in that 
case is not directly related to the black hole abundance, and the black hole bound on 
^(fcpeak) is unlikely to be strengthened very much. For instance, if the observational 
bound on black hole abundance applies separately to the black holes formed within 
each unit interval of In k, the effective value of y for a given value of (fc pe ak) is just 
multiplied by that factor, which has a negligible effect on the bound on V^(k peak ). 

4.2 The effect of Cm 

Now we discuss the effect of Ciim assuming that it is at least not canceled by Cni- By 
virtue of Eq. ( |2.5| ), the first term of Eq. ( |3.12| ) is <C 1, and the second term is < 1. 
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If /c*(t en d) is super-horizon, ( is of the form Eq. (|4.5|) with the minus sign, and the 
black hole bound is < (fc*(i e nd)) ^ 2. This is likely to be well satisfied. 

If instead fc*(i e nd) is sub-horizon, we have to remember that the black hole bound 
refers to horizon-sized regions. To apply it, we must drop sub-horizon modes of (jin- 
Estimating the bispectrum, trispectrum as in J| , one sees that this makes (j in nearly 
gaussian. Then P^ lin peaks at k ~ A; cn d = a(t CQ d)H, and the black hole bound is 
roughly P^ lin (/c e nd) < 1CT 2 . This too will be satisfied if fc*(t en d) is well within the 
horizon. 

We emphasize that these bounds refers to the formation of horizon-sized black 
holes. If fc*(t e nd) is sub-horizon, smaller black holes may also be formed. A discussion 
of their abundance would require assumptions about the evolution of the perturba- 
tions during the transition from inflation to non-inflation, and would be much more 
difficult than the corresponding discussion [R2J for the formation of black holes from 

Although V^ lin (k) is probably too small to form black holes, it may still be quite 
large. If reheating after inflation is long delayed this may lead to copious structure 
formation with a variety of possible cosmological effects [33|. 

Finally, let us see whether V^ lin (k) can be significant on cosmological scales; ie. 
whether it can be comparable with the observed quantity — 10~ 9 . It follows 
from Eq. ( j2.15 ) that the scale fc*(i e nd) is shorter than the scale leaving the horizon 
at the beginning of the waterfall. Therefore, the inequality ( p,12| ) implies that P^ lin 
will give a negligible contribution to the observed quantity ~ 10~ 9 if the shortest 
cosmological scale leaves the horizon more than 31n(10) ~ 7 e-folds before the start 
of the waterfall, ie. if the observable universe leaves the horizon more than ~ 22 
e-folds before the start of the waterfall. We will see that this is assured within the 
approximation of Section |[ 



5. Estimates using a simple approximation 

In this section we make a simple approximation for m 2 (t). This will allow us to verify 
some of the assumptions that we have been making, especially if we assume that <p 
satisfies the slow-roll approximation. 

5.1 The approximation for m 2 (t) 
The approximation is 

m 2 (t) ~ -/A (0 < iiH < m 2 ) (5.1) 
m 2 (t) ~ -m 2 (/A > m 2 ), (5.2) 
H 3 = 2#m|0(O)|. (5.3) 
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The cross-over between the two expressions is at t — t= = m 2 //i 3 . The second 
expressions corresponds to setting = 0. If the linear era ends at t < t= only the 
first approximation is invoked. 

The first approximation is exact at t — 0, and it ignores the time-dependence of 
d((p 2 )/dt = 200. The constancy of is a good approximation at t <C £=, and so is the 
constancy of if ( |5.8| ) is sufficiently well satisfied. The fast transition requirement 
described in Section [L3| is H < ji. To simplify some of the estimates we will usually 
take the requirement to be 

H <C fi (fast transition). (5-4) 

With this approximation for |m 2 (t)|, the linear era is completely described by 
the four parameters g, H, m, and fi. Let us define N(t) = Ht. Then the epoch 
t = t= corresponds to 

m\ 2 H ( m\ 3 H /m\ 2 [H x?J 



»M ={j)7= UJ » = W [7) ' (5 ' 5) 

5.2 Slow-roll approximation 

To obtain the strongest possible results, we assume that the evolution of satisfies 
the slow-roll approximation, at least during some era that begins before the waterfall 
and ends when ceases to to affect the evolution of x- 
Then unperturbed inflaton field satisfies 

0(t) + 3#OO0(t) + V'((f>(t)) = 0. (5.6) 

The basic slow-roll approximation is 

3#0~ -V ((/>), (5.7) 

or equivalently 

#|0/0| < 1. (5.8) 



The requirement that the first derivative of Eq. ( |5.7|) be consistent with Eq. (|5.8j ) 
is ej^ + ?7 <^ 1 where en = \H\/H 2 and r\ = V"/3H 2 . The slow-roll approximation 
assumes en «C 1 and \r)\ 1. 

Before the waterfall, and for as long afterward as x nas a negligible effect on 0, 
the Fourier components of the perturbation 50 satisfy 

50 k + 3#(t)50 k + (k/a) 2 54> k + K"(0(t))50 k = 0. (5.9) 

This equation ignores back-reaction, which is a good approximation by virtue of the 
slow-roll approximation H. 
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As a scale leaves the horizon, the vacuum fluctuation of <fi is converted to a 
classical perturbation 5<p with spectrum ~ (H/2tt) 2 . At a given epoch, the vacuum 
fluctuation is set to zero on sub-horizon scales. 

These results hold both before and during the waterfall. Focusing on the former 
era we have more results, because is the only field. First, we have a couple more 
relations: 

3M 2 H 2 (t) = P4> ~ V = V + V{4>) (5.10) 
^ = ^--Mi(^)72. (5.11) 

Second, we have the crucial result that 8(f) generates nearly gaussian curvature per- 
turbation to the curvature perturbation (, with spectrum given by 



H 2 



27T0 



(5.12) 

aH=k 



For a given k the spectrum is generated at the epoch of horizon exit k = aH, and is 
constant thereafter until at least the beginning of the waterfall. 

5.3 Trading \i for / 

For single- field inflation, we can use Eq. ( 5.12f) to obtain more powerful results by 



trading \i for |§] 

/ = (5 x lO- 5 ) _1 i? 2 /27T0(O) = (5 x lO- 5 )" 1 ^ 2 ^), (5.13) 

where fcbeg is the horizon scale at the beginning of the waterfall. Inflation models are 
usually constructed so that accounts for the observed on cosmological scales. 
Then, if is nearly scale-independent we will have / ~ 1. More generally there is 
an upper bound 

/ < 2 x 10 3 (black hole constraint) (5.14) 

corresponding to the black hole bound Vq, < 1CT 2 on the spectrum of the nearly 
gaussian C = C</> that exists at the beginning of the waterfall. There is also a lower 
bound corresponding to Eq. ( |1.6|) : 

/ > 10~ 2 #/(l(T 5 Mp) (nearly exponential inflation). (5.15) 

The relation between / and /i is given by 



— ~ HT 4 ^ — .. (5.16) 
\fij gm 

We are demanding g<l, but the fast transition requirement H \x can always be 
satisfied because / < 2 x 1CT 3 and H m. 
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In this paper we are not specifying the potential V(<f>). Most previous work 
considers the potential V(0) = |m 2 2 . Then slow-roll requires <C H and 3H<p = 
— m?0. The fast transition requirement H < fi becomes 

§)m>^ («■») 



and / is given by 

f = ^g[-) (-) . <r,it 



2 \ -1 



In this case we need / < 1, to avoid a positive spectral tilt for which would 
conflict with observation. The requirement that V(4>) does not support inflation (so 
that inflation ends with the end of the waterfall) is <p c < lOMp, which is guaranteed 

byEq. drrg). 

5.4 The case t end < t = 

In this case m 2 (t) ~ — fi 3 t, and Eq. ( |2.26[ ) becomes 



r 3g/i 3 # 2 

\ 2g 2 m? ' g ; 



Xni-min<{ ^f^,^ [>• (5.19) 



Using Eq. ( |2.26|) and x 2 ~ |?n 2 (£)|x 2 



-r— — ~ mm <^ 6iV(t end ),4iV(t end ) — 

2 tend I ^ 



(5.20) 

Growth begins when Eq. ( |2.6| ) is satisfied, corresponding to \xt\ ~ 1, and //i s tart = 
2y,tx ~ 2 and Ht staxt < 1. 

The case iV(t end ) «C 1 is considered in M. We then have 

kl(t end ) = a 2 (t start )^ 2 /2(M end ) 1/2 . (5.21) 
Using Eqs. ( gig ) and (p|) this gives* 10 Xni - /i 2 exp[(4/3)(/it end ) 3 / 2 ], ie. 

^ end ~ An ^ 7 . (5.22) 



Assuming instead iV(t end ) > 1, 

kl(t cnd ) ~ a 2 (* star tV 3/2 # 1/2 . (5.23) 
Using Eqs. (|2.18|) and (^8| ) we again arrive at Eq. (|5.22|) #11 



# 10 In the prefactor of Eq. ( |2.18 ) we drop a numerical factor and a factor t -3 / 2 , because these are 
negligible compared with the exponent. 

# n Factors fi/H are ignored in the prefactor of Eq. (2.18) because the exponent is > n/H. 
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Our implicit assumption that the growth era starts well before t en d corresponds 



to /itend 3> 1. This is equivalent to Xni 3> \x or 



mm 



H 



imgf 1 / 5 



5/3 



g 4 m l 



2/3' 



> 1. 



(5.24) 



The requirement t en d < £= corresponds to 

A^end < 



mN 2 
H 



10" 



m\ 2 f 



2/3 



(5.25) 



To get an upper bound on /it en d we use \i > H and xti ^ Xvev ^ to find 
^end ^ ( m (Mp/if)) 2 / 3 . Using the range ( |1.9[ ) this gives an upper bound on /it en d 
that is of order 5 to 20 with the lower end of the range far more likely. The upper 
bound corresponds to iV(t end ) < (H / /i)(ln(Mp / H)) 2 ^ 3 giving 



N(t cnA ) < [ln(M P /if)] 2 / 3 < ln(M P /F). 



(5.26) 



Using (x 2 (tst ar t)) ~ /i 4 and Eqs. {^J§, (^TJ, and |[24|),# 12 we find (x 2 (Wt)> < 
2 (^start)- This justifies the assumption |5i(x, tend) I ^ |<5t(x, t sta rt) | made after Eq. Q3.7| ) 
because the initial strong growth of (x 2 {t)) /<p 2 (t) will outweigh the slower variation 
of the other factors. 

5.5 The case t cn d > £= 



In this case m 2 (t en d) 



-m 



As we discussed in Section |3~1"| x 2 i is given by Eq. ( |2.25| 



and 7\ lin by Eq. (|3.15 ). Growth starts, at the latest, at t= + m -1 .* 13 

Our our approximation makes d\m(t)\/dt discontinuous at t — t= in violation of 
the adiabaticity condition (2.5). In reality |m(t) | will be smooth around t = t = . To 
avoid specifying a definite form for |m(i)|, we confine ourselves to the case N(t en d) > 
1. Then 

A; 2 (t end ) ~ a 2 (t start )m#. (5.27) 
Assume first that the growth era starts before t = . Then Eqs. (|2.18|) and 



#14 



give 



We also have 



H 2 exp ^2m(t cnd -t=) + ^ (jit = 



.3/2 



mt- 



(/it=) 2/3 = (m//i) 3 > 1, 



(5.28) 



(5.29) 



# 12 Only the first case of Eq. (5.24) need be invoked for this purpose. 

^ 13 Since growth always begins, \ always becomes classical. We have no need of a purely quantum 
treatment, which would require an entirely different approach. (In [Q we noted \ may fail to become 
classical within the regime t <C t—. Contrary to what was stated there, that does not imply that \ 
may fail to become classical at all.) 

# 14 Since mi enc j 3> m/H, and we ignore factors of m/H in the prefactor in Eq. (|2.1§|). 
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where the inequality holds because we are assuming that growth starts before t = 
Using the first equality we get 



XL ^ H 2 exp \2m U end - V)) J ~ H 2 exp (2mt cnd ) . (5.30) 

The final approximation is t= <C 3t en d, which should be adequate because we are in 
the regime t = < t end . In this approximation, the growth before t = has a negligible 
effect. Using it we find 

X 2 vev ~H 2 e 2mt ^, (5.31) 

leading to* 15 

N(t cnd ) = -(mt end ) ~ - ln(M P /#). (5.32) 
m m 

This gives again the bound ( |5.26| ). 

Now suppose that growth does not start before t = . Then the inequality in 
Eq. ( |5.29| ) is reversed leading to N= <C 1. We therefore arrive again at Eq. ( |5.31| ) 
leading to Eq. ( |5.32| ). In this case Ht start = N= + H/m which is <C 1 as before. Also, 
from Eqs. ( [2.1 3|) and ( |5.27|) , we have K(t staT t) / K(t end ) ~ H/m. Using Eq. ( j3T8[ ), this 



ensures that the typical value of |5£(x, t end )/5t(-x, t sta rt)| is > m/H 3> 1- 
5.6 Duration of the non-linear era 



For case t cnd > t = , xl nd is the value of x 2 at which the right hand side of Eq. ( |1.12j ) 



becomes important, corresponding to ln(x vcv /x C nd) ~ 1- For the opposite case, 
Eqs. ( |5T9D and (|5T|) give 

Wv <^. (5.33) 



v 2 M 2 

Avev ^ IV1 P 
Xend mH ' 



This is much less than Mp/H 2 , which means that Eq. (|3.16| ) gives iV n i ln(Mp/iJ). 
This is the same bound that we obtained for N(t end ). It therefore applies to the total 
number of e-folds of the waterfall, iV watcr ee iV(t end ) + iV nl . 



As seen in Section [4.2| , we need the waterfall to begin more than 22 e-folds after 
the observable universe leaves the horizon, if we are to be sure that V^ Un (k) has a 
negligible effect on cosmo logical scales. Equivalently we need N ohs — 22 > iV wat . From 
Eq. ( |1 . 1| ) the left hand side of this inequality is bigger than 47— [ln(Mp/if )]/2 and we 
have seen that the right hand side is -C ln(Mp / H) . The inequality will therefore hold 
if 47 > [ln(M P /#)]/2 ie. if H/M P > 10~ 41 . This is hardly stronger than the BBN 
bound ( |1.9|) , which means that V^ lin (k) is almost certainly negligible on cosmological 
scales. 

5.7 Two inflation models 

To illustrate the power of our results we apply them to two inflation models based 
on supersymmetry. 



# 15 We ignore a factor H/m within the log, which is permissible since H/m is also the prefactor. 
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5.7.1 Supersymmetric GUT hybrid inflation 



Supersymmetric GUT hybrid inflation |L2|, [L3] takes x to be a GUT Higgs field so that 
Xvcv — 10~ 2 Mp corresponding to (H/m) 2 ~ 10~ 5 . This is not small enough for the 
'end of inflation' formula to yield the entire waterfall contribution to ( (Eq. ( |3.22|) ). 
Supersymmetry gives g 2 = 2 A leading to g 2 = 10 9 (H/Mp) 2 . This leaves for our 
discussion two independent parameters which we take as g and /. 

The potential V((j)) may depend on several parameters |T3[. It typically steepens, 



and our discussion applies only if the parameters are such that steepening does 
not end slow- roll before t e nd- Requiring the inflaton perturbation to generate £ on 
cosmological scales, the steepening implies 10~ L5 g < / < 1, the lower bound coming 
from Eq. Q5.15|) . Using Eq. (|5.16|) we have [i/H ~ 10 2 (g/ f) 1 ^ 3 . The fast transition 



requirement (|5.16| ) is certainly satisfied if g ^> 10 (ie. H/Mp ^> 10 which 



usually holds [13]. 



The parameter space allows t cn d < t = (with either of the possibilities in Eq. ( |2.26|) ) 
as well as t en d > t = . Provided that Hj [i is well below 1, the duration of the waterfall 
is quite short, and the 'end of inflation' formula can give £ni in part of the parameter 
space ((ET2ID). 

5.7.2 Supernatural/running-mass inflation 



Supernatural inflation |14j and running-mass |L5| inflation take Xvev roughly of order 
Mp corresponding to m roughly of order H. This can be motivated by supposing that 
X is a string modulus, with gravity-mediated or anomaly-mediated supersymmetry 
breaking J34]. The former case, U 1/4 ~ 10 10 GeV orff~ 10~ 15 Mp is usually invoked 
and the latter would give H ~ 10~ 13 Mp or so. This low inflation scale and m ~ H 
are distinguishing features of the paradigm. 

The potential for supernatural inflation is V(<fi) = m|0 2 /2 which does not allow 
C = (</> on cosmological scales. Running-mass inflation takes V{<ft) to be the renor- 
malization group improved potential allowing ( = on cosmological scales which 
is assumed. In a suitable regime of parameter space, Cf>(k) on small scales can be 
big enough to exceed the cosmological bound on black hole formation, providing a 
constraint on the parameter space; in other words we can have / ~ 10 3 . This is 
another distinguishing feature of the paradigm. 

Since m is roughly of order H our criterion m 2 ^> H 2 cannot be very well 
satisfied and the analysis of the next Section is really more appropriate. To proceed 
we assume that m/H is a bit above 1, and take / ~ 1. Then the fast transition 
requirement H <C fi is satsfied for g 2 3> 10~ 8 which is expected. 

Since m is roughly of order H and we deal with the case t en d > t=, corresponding 
to m 2 (t en d) — — m 2 . This gives N(t en d) ~ ln(Mp/if) ~ 33, and fc*(t e nd) is outside 
the horizon, with V^ lin (k^(t en( j)) — (H /2m) 2 . This is is not far below 1, and the black 
hole bound might be violated. 



- 25 - 



The duration of the non-linear era is N n \ ~ 1. Since Eq. ( [3.21 ) is not satisfied, 
the contribution ( n \ is not given by the 'end of inflation' formula. 



6. The case m ~ H 

Now we consider the case that m/H is < 1 but not extremely small. To arrive 
at estimates we assume that V x {k) in this regime continues to peak at some value 
k*(t) <C a(t)\m(t)\. Since \m{t)\ < m ~ H this means that k*(t) is always outside 
the horizon. 

We assume that the gradient of x is negligible, checking the self-consistency 
of that assumption later. Then Eq. (|2.38|) holds. Considering either of the two 
independent solutions we define s(t) by 

X(*,t) = Hs(t) X (x,t), (6.1) 

giving 

s 2 (t) + 3s(t) - \m 2 {t)\/H 2 = -s/H. (6.2) 



We assume that the right hand side of Eq. (|6.2|) is negligible, checking later for 
self-consistency. Keeping only the growing mode this gives 



3 9 \m 2 (t)\ . s 

^ = -2 + V4 + 4^ (0) 



l\m 2 (t)\ 
3 H 2 



(6.4) 



The gradient of x is indeed negligible compared with x, which means that p x and p x 
are given by Eqs. (|2.29|) and ( |2.31| ). Also, the assumption that the right hand side 



of Eq. ( |6.3| ) is negligible is self-consistent if and only if Eq. ( |2.40| ) holds. This means 
that the field equation Q2.38| ) is consistent with the energy continuity equation. 

If the approximation ( |5.1|) holds, Eq. ( p.40|) becomes Ht ^> 1. But the same 



approximation inserted into Eq. (|6.1|) gives at Ht ^> 1 

Ht ~ {H/pf In [ X 2 (x, t)/ X 2 (x, H- 1 )] . (6.5) 

Since we haven't calculated x( x 5^) from the vacuum fluctuation we don't know the 
precise value of x 2 ( x > H~ l ) but it presumably lies roughly between \i and m ~ H 
since these are the relevant mass scales. As we saw earlier this would make the log 
at most of order 10 2 or so. Therefore, since we are imposing H /i, Eq. ( |6.5|) is 
hardly compatible with Ht ^> 1. The conclusion is that Eq. (|2.40| ) probably requires 
the regime ( |5.2| ), m 2 (t) ~ — m 2 , which we assume from now on. That in turn implies 

(X 2 (4nd)) > <j> 2 (t cnd ). 
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To calculate £iin we use Eq. (3.4), and assume \5t(t CQI i) / 5t(t star t)\ 1- Using 
Eqs. ( ggp , (pD, and Q, we find 

p x /f = 3/2s(t) ~ # 2 /™ 2 , (6-6) 
to be compared with p x /x 2 = 3H/2\m(t)\ in the case m 2 3> i? 2 . We therefore have 



Clin(x) 



H 2 5 X 2 (x,tend) 



(6.7) 



m 2 (X 2 (tend)} 

Since V x (k) peaks at fc*, we expect that the final equality of Eq. fl2.18|) will be 
roughly correct. Also we expect that Vs x a{k) will be given roughly by Eq. (|2.22|) at 
k < fc* and will fall off at bigger k. Then, using Eq. ( |6.4| ) we see that V^ lin {k) peaks 
at k ~ &*(t e nd) with a value 



7VO*0end)) 



# 2 



(6i 



We conclude that the black hole bound is likely to be violated if m is significantly 
below H. 

A crucial feature of our setup is the condition ( [2.40 ), which is necessary for 
consistency if the gradient of x is negligible and there is no cancellation between the 
two terms of p x . We found that the solution of Eq. ( |2.1| ) then indeed makes the 
gradient of x negligible, with no cancellation. But the solution of Eq. (2~1) may also 
make the gradient negligible with no cancellation, in a part of parameter space where 
with the condition ( |2.40| ) violated. In such a regime, we would have to conclude that 
the linear approximation leading to Eq. (|2.1| ) is invalid. 



7. Comparison with other calculations 



Nineteen other papers have considered the contribution of the waterfall to x [14, 35 



36], |37|, [38], |9[[|y, H, (H ||@ H, p, ||. 11011], |5|, H . in the fast transition 



regime. Some of them also consider the issue of black hole formation fi"4] , 35 , |5T] , |53] , 
concluding that the black hole constraint is satisfied for m ^> H but not for m ~ H . 
That is roughly our conclusion though we are less sure. In view of this, one may 
wonder whether the present paper and its companion M are needed. They are, for 
several reasons. 

First, all of the previous papers take <fi to be canonically normalized, and nearly 
all of them go much further by assuming V(4>) = m 2 </> 2 /2. Second, none of the 
previous papers specifies all of the assumptions that are made, as we do here and in 



[[| . Third, none of them except [^1], [53|] considers the non-gaussian black hole bound 
as we have done in the present paper. Fourth, most of them present a calculation 
which is much more complicated than ours. Finally, all of the other papers except 
perhaps [[HI have errors. 
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The last point was considered in our earlier paper H. The problem for many 



of the papers [14], |35], ^6], [37], |8|, |39], [40], [41], 42] is that the waterfall is treated as 



two-field inflation, without imposing the requirement (x 2 ) (%>) that would be 
needed to justify such a treatment.* 16 As we have seen, this is not the case. It is 
only within the slow transition regime, considered in [52], [54], [55], [56|], that one can 
expect to find a regime of parameter space that allows the waterfall to be treated as 
two-field inflation. 



We will not repeat the analysis of the problems of the other papers IfiJ [45], [46 
[49] , |50ll appearing before |4j], that was given in the latter paper. After 



47| , [48| , [49| , |50l appearing before |4J, that was given in the latter paper. After [|4| three 
papers appeared [[40], |5T, [53|. The paper [[4(J is a continuation of [45]. Its main focus 
is on the case (f) c > lOMp, where inflation continues after the waterfall, but that does 
not affect the contribution to ( x generated during the waterfall, and Eq. (6.15)) of 
[|4Cfl reproduces the expression for ( x given in Eq. A(10) of fl4"5"f . 



The papers [^1], consider the case m < H, and they calculate \k by numer- 
ical integration with the potential V(<f)) = m^0 2 /2.# 17 Then they evaluate (j in by 
integrating Eq. fl3.2| ), finding that m < H is definitely forbidden by the black hole 
bound. The calculation assumes that the gradient of % can be ignored when evalua- 
tion p x and p x but they don't investigate the compatibility of the evolution equation 
Eq. (|2.1|) with the energy continuity equation. However, their results for the case 
m = H (with the other parameters fixed at particular values) shown in their Figure 
2 is in excellent agreement with our Eq. ( |6.1|) , assuming |m(t)| = m. Their result for 
^ > Ciin(^*(^end), ^cnd) with the same parameter choice, shown in their Figure 4, is also in 
agreement with ours, assuming in addition (x 2 (£ en d)) 3> 2 (^end)- It therefore seems 
that for at least this parameter choice, their assumption that the gradient of x is 
negligible is justified, and that moreover the consistency condition (|2.40|) is satisfied. 
But there is no reason to think that the same is true in the entire parameter space, 
considered in their Figure 5. Regarding the black hole bound, they note that (n n 
has the non-gaussian form ( |4.5|) . In |5l| they use (£ 2 n ) < 1 instead of our V^ hn < 1. 
As the width of the peak in V^ Un (k) is rather broad, this will somewhat overestimate 
the region of parameter space forbidden by the cosmological bound on V^ lin as we 
noted in Section (4~T[ . In |53[ a more sophisticated procedure is used to obtain the 
black hole bound, but they don't estimate the theoretical error and it is unclear to 
us whether it represents an improvement on our rough estimate < 1. 



8. Conclusion 

We have considered the contribution £ lin to £, that is generated during the linear era 
of the waterfall within the Standard Scenario. We gave a rather complete calculation 



# 16 The papers (3^, 39, set \ 2 = {\ 2 ) while the others regard it as a free parameter. 



# 17 With m < H, the fast transition requirement (5.17) conflicts with the slow-roll requirement 
Tfifp <§; H, but the two are roughly compatible with the choice (m^/H) 2 = 10 _1 of |5l], p^ . 
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for the case that the waterfall mass m is much bigger than if, and arrived at estimates 
for m ~ H . 

Taking on board our discussion of the non-gaussian black hole bound, we con- 
cluded that the black hole bound will be satisfied for m ^> H, but that it may well 
be violated for m < H . The latter case will be further investigated in a future 
publication [[5711 , by numerically integrating Eq. Q2.1|) . 

A lot more will have to be done before we have a complete understanding of 
contribution to ( generated during the waterfall. A fundamental problem is to handle 
the ultra-violet cutoff, that is needed to obtain finite values for the fields and for the 
energy density and pressure. Our procedure of keeping only the classical field modes 
is approximate, and it violates at some level the energy continuity equation. This 
and related issues are discussed for instance in [0 . A precise procedure is advocated 
in p9| , but its relation to our approximate procedure is unclear. 



An understanding of the ultra-violet cutoff will allow one to decide on the min- 
imum value of 0(0) that allows an initial linear era [17| , |J. With that in place one 
would hopefully verify that the value invoked in the present calculation is big enough. 
But it will still be unclear how to evaluate the contribution to ( that is generated 
during inflation after the linear era ends, when it is not given by the 'end of in- 
flation' contribution. A numerical simulation, even with reasonable simplifications, 
might well require one to consider a patch of the universe that is too big to handle. 
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